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Abstract 

>, 

^P ■ We present a measurement of the CKM matrix element \Vch\ and the form- factor slope p^ for 

Q>^ \ B -^ Di~i'£ decays based on 417 fb~^ of data collected at the T{AS) resonance with the BABAR 

■^ ' detector. The semileptonic decays are selected in BB events in which the hadronic decay of the 

t~^ ■ second B meson is fully reconstructed. From the measured differential decay rate of the signal 

O ; decay we determine g{l)\Vcb\ = (43.0 ± 1.9 ± 1.4) x 10"^, p'^ = 1.20 ± 0.09 ± 0.04, where ^(1) is 

^— s I the hadronic form factor at the point of zero recoil. Using a lattice calculation for ^(1) we extract 

\ycb\ = (39.8 ± 1.8 ± 1.3 ± 0.9) X 10^^, where the stated errors refer to the statistical, systematic, 
and form factor uncertainties. We also present a measurement of the exclusive branching fractions, 
X ; B{B- -^ D'^i-ve) = (2.31 ± 0.08 ± 0.07)% and ^(5° -^ D+rvi) = (2.23 ± 0.11 ± 0.( 
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1 Introduction 

The study of the semileptonic decay B -^ Di~v'£ is interesting in many aspects. In the standard 
model, the rate of this weak decay is proportional to the square of the Cabibbo-Kobayashi-Maskawa 
(CKM) pj matrix element \Vcb\, which is a measure of the weak coupling of the b to the c quark. 

The decay rate is also proportional to the square of the hadronic matrix element, which accounts 
for the effects of strong interactions in the B ^ D transition. In the limit of very small masses of 
the lepton, i = e or fi, their effect can be parameterized by a single form factor G^w), where the 
variable w (see below) is linearly related to the momentum transfer squared q^ of the B meson to 
the D meson. 

The extraction of \Vcb\ relies on the measurement of differential decay rates of semileptonic B 
decays. The precise determination requires corrections to the prediction for the normalization at 
w = 1 in the context of Heavy Quark Symmetry [2j , as well as a measurement of the variation of 
the form factors near the kinematic limit, w = 1, where the decay rate goes to zero as the phase 
space vanishes. 

Measurements of \Vcb\ based on studies of the differential decay rate for B -^ Di~i>i decays 
have previously been reported by Belle [3], CLEO [4j and ALEPH |5j. 

In this paper, we present preliminary measurements of the differential decay rates, separately 
for B^ -^ D~^£~r'£ and B^ -^ D^l^vi [6]. Semileptonic decays are selected in BB events in 
which an hadronic decay of the second B meson is also fully reconstructed. This leads to a very 
clean sample of events and also provides a precise measurement for the variable cf' , and thereby w. 
The measurement of the total branching fractions and the extraction of \Vcb\ requires an absolute 
normalization. We use a sample of inclusive semileptonic decays, B -^ Xi'v'i, where only the 
charged lepton is reconstructed. This choice reduces the uncertainty of the normalization, because 
the lepton is also selected from a sample tagged by hadronic decays of the second B meson. 

2 Parameterization of the Decay Rate 

In the limit of small lepton masses the partial decay rate for B -^ Dl~U£ can be expressed in terms 
of a single form factor, 0(10), 

where Gp is the Fermi coupling constant, and Mb and Af/j are the masses of the B and D mesons. 
The variable w denotes the product of the B and D meson four- velocities Vb and Vd, 

{Ml + Ml- q') 
^-^b-Vd- ^2MbMo) ' ^^^ 

where q^ = {ps — Pd)'^-, and pb and pn refer to the four-momenta of the B and D mesons. Its 
lower limit, w = 1, corresponds to zero recoil of the D meson, i.e. the maximum q^. The upper 
limit, w = 1.59, corresponds to g^ = and the maximum D momentum. Since the B momentum 
is known from the fully reconstructed -Btag in the same event, w can be reconstructed with good 
precision, namely to ~ 0.01, which corresponds to about 2% of the full kinematic range. 

In the limit of infinite quark masses Q{w) coincides with the Isgur-Wise function [^. This 
function is normalized to unity at zero recoil. Corrections to this prediction have recently been 
calculated with improved precision, based on unquenched lattice QCD [8], specifically ^(1) = 



1.074 lb 0.018 lb 0.016. Thus \Vcb\ can be extracted by extrapolating the differential decay rates to 
w = 1. To reduce the uncertainties associated with this extrapolation, constraints on the shape of 
the form factors are highly desirable. Several functional forms have been proposed [9|. We adopt 
the parameterization based on analyticity and positivity of the QCD functions which describe the 
local currents [1^. Specifically, G{w) is expressed as a polynomial in z, 

g{w) = g{l) h - Sp'^z + (51p2 - 10)^2 - (252p2 - 8i)z^] , (3) 

where the variable z is defined as z = {\/w + 1 — v2)/(\/uj + 1 + v2)- This formulation expresses 
the non-linear dependence of the form factor on w in terms of a single shape parameter, p^. 

3 The BaBAR Detector and Dataset 

This analysis is based on 417 fb^^ of data collected at the T(45) resonance with the BABAR detector 
at the PEP-II storage rings. The corresponding number of produced BB pairs is 460 million. In 
addition, 40 fb~^ of data, recorded at a center-of-mass energy 40 MeV below the T[4S) resonance, 
are used to study background from e~^e~ -^ ff {f = u,d,s,c,T) events (continuum production). 
The BABAR detector is described in detail elsewhere |llj . Charged-particle trajectories are measured 
by a 5-layer double-sided silicon vertex tracker and a 40-layer drift chamber, both operating in 
a 1.5-T magnetic field. Charged-particle identification is provided by the average energy loss 
{dE/dx) in the tracking devices and by an internally refiecting ring-imaging Cherenkov detector. 
Photons are detected by a CsI(Tl) electromagnetic calorimeter (EMC). Muons are identified by the 
instrumented magnetic-flux return. A detailed GEANT4-based Monte Carlo (MC) simulation |12j 
of BB and continuum events has been used to study the detector response, its acceptance, and 
to test the analysis procedure. The simulation models the signal B -^ Di~U£ decays using the 
ISGW2 model [13], and these are then reweighted to the HQET model described above. Other 
semileptonic decays that contribute to background are simulated as follows: for B — > D*(.~i'i 
decays we use form factor parametrizations based on Heavy Quark Effective Theory (HQET) with 
parameters determined by the BABAR collaboration [1^, for B -^ D**{-^ D^*>7r)£~i'g decays we 
use the ISGW2 model [13j, and for decay involving non-resonant charm states, B — > D^*'7r£~i'£, 
we adopt the prescription of the Goity-Roberts model [15] . The MC simulation includes radiative 
effects such as bremsstrahlung in the detector material. QED final state radiation is modeled by 
PHOTOS [TBj, and decays with radiative photons are included in the signal sample. 

4 Event Selection 

We select semileptonic B meson decays in events containing a fully reconstructed B meson (Stag), 
which allows us to constrain the kinematics, to reduce the combinatorial background, and to de- 
termine the charge and flavor of the signal B. 

The analysis exploits the presence of two charmed mesons in the final state: one is used for the 
exclusive reconstruction of the i?tag! and the other one for the reconstruction of the semileptonic 
B decay. 

We first reconstruct the semileptonic B decay, selecting a lepton with momentum in the center- 
of-mass (CM) frame p^ higher than 0.6 GeV/c. Electrons from photon conversion and vr^ Dalitz 
decays are removed using a dedicated algorithm, which performs the reconstruction of vertices 
between tracks of opposite charge whose invariant mass is compatible with a photon conversion or 



a n^ Dalitz decay. Candidate D^ mesons, with the correct charge correlation with the lepton, are 
reconstructed in the K~7r~^, K~tt^'7T^, K~ t:~^ tt'^ t:~ , K^t[^t[~ , K^'k^'k~'k^, K^-k^ , K~^K~, tt+tt", 
and K'^K'^ channels, and D+ mesons in the K~7r+TT+, K~TT+ir+TT^, K'^tt^ , A'^7r+7r°, K+K-ir+, 
K^K~^, and K^7r~^7r^TT^ channels. D candidates are selected within 2a of the D mass, with a 
typically around 8 MeV/c^. In events with multiple B -^ D(.~Di candidates, the candidate with 
the largest D-i~ vertex fit probability is selected. 

We reconstruct -Btag decays of the type B -^ D^*'Y, where Y represents a collection of hadrons 
with a total charge of ±1, composed of nivr + n2K + n-^Kg + n^ir^, where ni + n2 < 5, n^ < 2, 
and 714 < 2. Using D^{D^) and D*^{D*^) as seeds for B~{B^) decays, we reconstruct about 1000 
decay modes. 

The kinematic consistency of a -Btag candidate with a B meson is checked using two variables: 

the beam-energy substituted mass thes = y ^/^ ~ Pbi ^'^'^ ^^^ energy difference IS.E = Eb — y/s/2. 
Here y/s refers to the total CM energy, and pB and Eb denote the momentum and energy of the 
Stag candidate in the CM frame. For correctly identified i?tag decays, the ruEs distribution peaks at 
the B meson mass, while /S.E is consistent with zero. We select a -Btag candidate in the signal region 
defined as 5.27 GeV/c^ < niEs < 5.29 GeV/c^, excluding i?tag candidates with daughter particles 
in common with the charm meson or with the lepton from the semileptonic B decay. In the case of 
multiple -Stag candidates, we select the one with the smallest |Ai?| value. For charged B events, the 
Stag and the Di.~ candidates are required to have the correct charge-flavor correlation. We do not 
apply any correction to account for B^ — B^ mixing effects, because they are found to be negligible. 
Cross-feed effects, i.e. B^ {B^^^ candidates erroneously reconstructed as a neutral (charged) B, 
are subtracted using the MC simulation. 

Semileptonic B decays are identified by the missing mass squared, defined as: 



m 



2 



miss 



b(r(4S))-p(Stag)-pp)-pW]' (4) 



in terms of the particle four-momenta. For correctly reconstructed signal events, the only missing 
particle is the neutrino, and the rn^^^^ peaks at zero. Other semileptonic B decays, like B — > D*£~Ui 
and B — > D**l~V(^, where one particle is not reconstructed (feed-down), spread to higher values of 
m^jgg. The use of the full-reconstruction technique results in a m^jgg resolution of 0.04 GeV^/c^, 
an order of magnitude better than achieved in non-tagged analyses |17] . 

5 Measurement of \Vcb\ and p^ 

We measure \Vcb\ and the form- factor slope p^ by a fit to the w distribution for B -^ Di~ve decays. 

Data and Monte Carlo events are collected in ten equal-size w bins. The few events that have 
w < 1.0 or It; > 1.59 due to resolution effects, are collected in the first bin and the last bin, 
respectively. 

To obtain the semileptonic B -^ Dd.^'vi signal yield in the different w intervals, we perform 
a one-dimensional extended binned maximum likelihood fit to the wi^jgg distributions, based on 
a method developed by R. Barlow and C. Beeston [18j. The fitted data samples are assumed to 
contain four different types of events: 

• signal B -^ D£~i>£, 

• feed-down semileptonic B decays. 
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• BB and continuum background, 

• fake lepton events. 

We use the Monte Carlo predictions for the different semileptonic B decay "^^iss distributions 
to obtain the Probability Density Functions (PDFs) to fit the data distributions. The feed-down 
contributions are allowed to float independently in the different w intervals. We use the off-peak 
data to provide the continuum background normalization, while the BB normalization is taken from 
the MC. The shape of the continuum background predicted by the MC simulation is consistent with 
the one obtained from the off-peak data. 

The JTi^jgg distributions for two different w intervals are compared with the results of the fits 
in Fig. [H 
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Figure 1: Fit to the w-^iss distribution, in two different w intervals, for B^ -^ D^£^i?£ (top row), 
B^ — >■ D~^(.~i>i (bottom row): the data (points with error bars) are compared to the results of the 
overall fit (sum of the solid histograms). The PDFs for the different fit components are stacked 
and shown in different colors. 

The signal yields in bins of w are displayed in Fig.[2j A y^ fit is performed on the w distributions 
comparing the number of events observed in each interval to the MC signal prediction. The signal 
prediction is obtained by properly weighting each MC event. The x^ is defined as: 
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-1 Kata)^ + Eft'^ W^f 

where the index i runs over the ten w bins; -/V^ata ^^ ^^^ observed number of signal events found in 
the i^^ bin and the cr^^^^^^ the corresponding error. The expected signal yields are calculated at each 
step of the minimization from the reweighted sum of iV^c simulated events. The Monte Carlo is 
corrected event by event for all known differences in tracking, cluster reconstruction and particle 
identification determined from data control samples. The weight for each Monte Carlo event is 
computed as the product of two terms W- = W'-' x W^' ^° where 

• W'-' is an overall fixed scale factor that accounts for the relative luminosity of data and signal 
Monte Carlo events; this is obtained, as described below, by normalizing to the inclusive yield 
of -B ^ Xi~v'i events, corrected for its reconstruction efficiency, and scaled to the absolute 
B — > X£^i>c branching fraction |19| . 

. W^'theo = f^--{wy,p'',g{l)\Va\)/fMc{wy,pl,c,Q{l)\Vcb\Mc) is the term that describes the 
differential decay rate. It depends on ^(l)|V"cfe| and p'^, and varies at each step of the mini- 
mization. 

Here w^, is the true w of the MC event j, reconstructed in the z-bin. The function f^^'^° corresponds 
to the expressions in Eq. [Hand El while /mc is the function used to generate B — > Di~r'i Monte 
Carlo events |j. 

We first fit separately the w distributions for the charged and neutral B -^ Di~v'i samples; we 

also perform a combined fit for B — > D£~i'i, assuming that /qo + /h = 1, i-e. that B^B~ and 

B^B^ together saturate the T(4S) decays. The value of the branching ratio is then computed by 
integrating the differential expression in Eq. [TJ 

In Fig. [2] we show the comparison between the data and the fit results separately for B~ — > 
D^i~i>l and B^ — > D^i~i'£ . The corresponding distributions for the combined fit are shown in 
Fig. [3l The measured values of ^(l)!^^^! and p^, with the corresponding correlation pcorr obtained 
from the fit, are reported in Table [Hand shown in Fig. [H 

For the measurement of ^(l)|T4b|, in order to reduce the systematic uncertainty, we use as 
normalization a sample of inclusive B -^ X(,~i>i decays. These events are selected by identifying 
a charged lepton with CM momentum greater than 0.6 GeV/c. In the case of multiple -Btag 
candidates, we select the one reconstructed in the decay channel with the highest purity, defined 
as the fraction of signal events in the tues signal region. We require the lepton and the -Stag 
to have the correct charge correlation and the lepton track to not be used to reconstruct the 
-Btag candidate. Background components peaking in the ruEs signal region include cascade B 
meson decays (i.e. the lepton does not come directly from the B) and hadronic decays where 
one of the hadrons is misidentified as a lepton. These backgrounds are subtracted by using the 
corresponding simulated Monte Carlo distributions. The cascade-i? meson decays (17.6% and 
19.0% of the total niES distribution for charged and neutral S, respectively) are reweighted to 
account for differences between the branching fractions used in our Monte Carlo simulation and 
the latest experimental measurements |20] . The total yield for the inclusive B -^ Xl~Di decays 
is obtained from a maximum- likelihood fit to the niEs distribution of the Stag candidates, using 



^Corresponding to the CLN parameterization pXij, with p^ = 1.17. 

12 



Table 1: Fit results for each sample. In the last column we report the results for the B^ and B~ 
combined fit, where the branching fraction refers to B^ decays. We also report the signal yields 
and the reconstruction efficiencies, integrated over the full w range. Only the statistical errors are 
reported here. 
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g{i)\Vcb\ ■ 10^ 
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Pcorr 

xVndf 
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Figure 2: The w distribution for B -^ D£ ug events. Left: B -^ D^£ ui, right: B^ -^ D^£ vi. 
The data (points with error bars) are compared to the results of the overall fit (solid histogram). 



an ARGUS function [21] for the description of the combinatorial BB and continuum background, 

and a Crystal Ball function [22] for the signal. A broad peaking component is observed in the ruES 
signal region and is included in the signal definition. This is due to real B — > Xi!.~h'i decays for 
which, in the -Btag reconstruction, neutral particles are not reconstructed or are interchanged with 
the semileptonic decays (e.g. a 7 from radiative D*^ decay which belongs to the D*^ seed in the 
-Btag decay chain and is instead associated with a B~ -^ D*^{D^'y)i~i'£ decay). This broad peaking 
component is modeled with additional Crystal Ball and ARGUS functions, whose parameters are 
fixed to the Monte Carlo prediction, except for the Crystal Ball mean value. Fig. [5] shows the 
ruES distribution for the i?tag candidates in the B~ -^ Xi^ui and B^ -^ Xi~^ij£ samples. The 
fit yields 198,897 ± 1,578 signal events in the B~ —>■ Xl~V£ sample and 116,330 ± 1,088 signal 
events in the B^ -^ Xi-~v^ sample. The corresponding reconstruction efficiencies, including the 
-Btag reconstruction, are 0.39% and 0.25%, respectively. 
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Figure 3: Left: w distribution obtained summing together B~ — > D^i~r'£ and B^ -^ D^(,~vi 
yields. The data (points with error bars) are compared to the results of the overall fit (solid 
histogram). Right: ^(1^)1^461 distribution unfolded for the reconstruction efficiency, with the fit 
result superimposed. These plots are not corrected for the smearing in w. 
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Figure 4: (left) Ax^ = 1 ellipses in the ^(l)|T4f,| versus p^ plane for the B^ -^ D^l^v^ decay, 
B^ -^ D^i.~V(^ decay and the B -^ Dl~vi combined fit (ellipse in red). The systematic uncertainties 
are taken into account. 



6 Systematic Uncertainties 

Different sources of systematic uncertainties have been estimated and are given in Table [2j We 
have grouped them into several categories. 

Detector-related systematics may arise from differences between the data and simulation of 
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Figure 5: niEs distributions of the a) B^ -^ Xi~V£, and b) B^ -^ Xi~vi samples. The data 
(points with error bars) are compared to the result of the fit (solid line). The dashed lines show 
the broad-peaking component and the sum of the combinatorial and continuum background. 

the track reconstruction and efficiency, particle identification and neutral particle reconstruction. 
The systematic uncertainty related to the reconstruction of charged tracks is determined by ran- 
domly removing a fraction of tracks corresponding to the uncertainty in the track finding efficiency, 
estimated on e^e~ — > t~^t~ data control samples. The systematic uncertainty due to the re- 
construction of neutral particles in the EMC is studied by varying the resolution and efficiency 
to match those found in data control samples. We estimate the systematic uncertainty due to 
particle identification by varying the electron and muon identification efficiencies by 2% and 3%, 
respectively. The misidentification probabilities are varied by 15% for both electrons and muons. 

We evaluate the systematic uncertainties associated with the Monte Carlo simulation of the 
various signal and background processes. The uncertainty arising from radiative corrections is 
studied by comparing the results using PHOTOS with those obtained with PHOTOS turned off. 
We take 30% of the difference as a conservative systematic uncertainty. The fraction of B cascade 
decays in the B -^ X£~iyi sample is varied within its uncertainties and the differences in the B — ^ 
Xi^i'i signal yields are included in the systematic uncertainties. Possible differences in the -Stag 
composition of the MC simulation and data can affect the efficiencies and the cross-feed between 
charged and neutral B events. To evaluate this effect we assume a conservative 30% systematic 
uncertainty to the cross-feed fractions and we evaluate the systematic uncertainty by looking at 
differences in the measured values of ^(l)|V'cb| and p'^ as we change the cross- feed fractions. We 
vary the B -^ D*i~iyi form factors within their measured uncertainties [U] and we use an HQET 



parameterization [23] for B — > D**i~v'£. We also vary the B — > D**(.~i>i branching fractions 
within the measured uncertainties |19j . For the \Vch\ measurement, we include a contribution due 
to the uncertainties on the branching fractions of the reconstructed D modes, and on the absolute 
branching fraction B{B — >■ X(.~i'() used for the normalization. 

We also evaluate a systematic uncertainty due to differences in the efficiency of the -Btag selection 
in the exclusive selection oi B ^ Di~ui decays and the inclusive B -^ Xi~v'i reconstruction, by 
using the same -Btag candidate selection adopted in the B -^ Xi^ug reconstruction also for the 
B — > Di~r'i decays, and taking the difference in the signal yield, corrected for the difference in the 
reconstruction efficiency, as a systematic uncertainty. 

The systematic uncertainty in the determination of the B -^ Xi.~i>i yield is estimated by using 
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Table 2: Systematic uncertainties in the measurement of t/(l) T^^ and / 


92 for B 


-^ Di^vi 


decays. 


We report the relative error (in %) for ^(1) 


Vcb and the absolute error on p'^. 








Systematic uncertainty on \Vcb 


and p^ 




dH- 


V£ 


D+t 


'Vi 


D£- 


i^e 




\Vc,\{%) 


p' 


\Vc,\{%) 


p' 


\Vch\{%) 


P' 


Tracking efficiency 


0.5 


0.008 


1.1 


0.003 


0.7 


0.004 


Neutral reconstruction 


1. 


0.003 


0.8 


0.006 


0.9 


0.004 


Lepton ID 


1.0 


0.009 


0.9 


0.009 


0.95 


0.009 


PHOTOS 


0.13 


0.005 


0.10 


0.005 


0.12 


0.005 


Cascade B ^ X ^ l'^ decay background 


0.6 


- 


1.0 


- 


0.75 


- 


B — B~ cross-feed 


0.24 


0.003 


0.24 


0.003 


0.24 


0.003 


B -^ D*£~i?e Form factors 


0.56 


0.008 


0.20 


0.003 


0.38 


0.006 


B -^ D**e~Ui Form factors 


0.24 


0.007 


0.34 


0.006 


0.29 


0.007 


D branching fractions 


1.0 


- 


1.35 


- 


1.12 


- 


B{B -^ D**i-ue) 


1.18 


0.023 


0.96 


0.011 


1.08 


0.019 


B{B -^ Xi~i)e) 


0.95 


- 


0.95 


- 


0.85 


- 


Bta,g selection 


1.1 


0.021 


1.8 


0.036 


1.5 


0.028 


B -^ Xi-Di yield 


0.7 


- 


1.1 


- 


0.85 


- 


B -^ Di~ue yield 


1.27 


0.018 


1.06 


0.027 


1.25 


0.020 


Total systematic error 


3.1 


0.04 


3.6 


0.05 


3.3 


0.04 



an alternative fit method, which is then compared to the result of the nominal niEs fit- We consider 
the niES distribution from the data and the combinatorial BB continuum and other background 
components (cascade and hadronic B decays) modeled with distributions taken from the Monte 
Carlo simulation. We fit the background normalization on data in the rriES sideband region, defined 
by rnES < 5.265 GeV/cP. The normalization for the continuum background is fixed to the value 
obtained from off-peak data. The total background contribution is then subtracted from the total 
number of events in the niEs distribution to extract the B -^ Xi^ui signal yield. The uncertainty 
in the determination of the B -^ Di^i/£ yield in the different w intervals is estimated by changing 
the PDFs used to model the different contributions in the m^^^^ distribution, e.g. by replacing the 
continuum PDFs with the corresponding one obtained from off-peak data. 

7 Results 

We present a measurement of ^(1)|V^5| and p^ for B -^ Di^v^. For i?~ -^ D^l~v^, we obtain: 



B{B- 



Qinvc 



cb\ 



P 



DH-Vi) 



(41.7 ±2.1 ±1.3) X 10-3 
1.14 ±0.11 ±0.04 
(2.31 ±0.08 ±0.07)%, 



(6) 



while for B 



■0 



D+t 



vi, 



we obtain: 
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^(1)1 


Vcb\ 




p' 


D+t 


'vt) 



g{l)\Vcb\ = (45.6 ± 3.3 ± 1.6) X 10"^ 
p2 = 1.29 ±0.14 ±0.05 
B{B^ ^ D+Tui) = (2.23 ±0.11 ±0.08)%. (7) 

The results of the combined fit are: 

: (43.0 ± 1.9 ± 1.4) X 10"^ 
: 1.20 ± 0.09 ± 0.04 
B{B^ ^D+rui) = (2.17 ±0.06 ±0.07)%. (8) 

Using an unquenched lattice calculation [5], corrected by a factor of 1.007 for QED effects, we get 

\Vcb\ = (39.8±1.8±1.3±0.9ff) X 10"^ (9) 

where the third error is due to the theoretical uncertainty in ^(1). The resulting value of \Vcb\ is fully 
compatible with the other existing measurements on i? — > Di~r'i, and also with the measurement 
obtained using B -^ D*i~i>i. Within the total errors this measurement is also compatible with the 
inclusive determination of \Vcb\ = (41.68 ± 0.39 ± 0.58) x 10"^ [24J. 
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